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We present a theory of the magnetic rare-earth nickel 
boride-carbides (RENi 2 B 2 C, RE=Ho, Er, Tm). Large total 
angular momentum allows one to employ semiclassical ap- 
proximation and reduce the problem to the four-positional 
clock model. The latter gives a proper phase diagram in mag- 
netic field-magnetic moment plane for varying directions of 
in-plane magnetic field. The theory explains recent experi- 
mental observations. 

PACS numbers 74.72.Ny, 75.30.Cr, 75.30. Kz, 75.10.Dg 



A complex phase diagram of HoN^I^C in the tempe- 
rature-magnetic field plane has been discovered in the 
neutron scattering experiments [jjj and in measurements 
of resistivity . Recent measurements of the anisotropic 
magnetization of H0NI2B2C by P. Canfield and cowork- 
ers ||,^| revealed a complicated magnetization-magnetic 
field diagram: re-entrant behavior for the low fields and 
a set of magnetic transitions for the high fields. The 
magnetic moment is highly anisotropic and has a strong 
dependence upon direction of magnetic field in the basal 
a-b plane. A purpose of this article is to present a simple 
model for the magnetic subsystem in RENI2B2C, explain- 
ing the experimental facts. 

The crystal structure of the RENI2B2C has been stud- 
ied by T. Siegrist et al. || (RE=Lu) and Q. Huang et 
al. § (RE=Ho). The RENi 2 B 2 C compounds have the 
body-centered-tetragonal (bet) structure (space group 
IA/mmm). A simplest form of the single- ion Hamilto- 
nian compatible with the crystal field symmetry is 



and J y = J sin <f>. Then Hamiltonian (Q) can be presented 
in the form 



Ucef = -^ J z + + J -)^ 



(1) 



where J±, J z are components of the total angular mo- 
mentum, a and b are crystal field constants. We assume 
that moments are presumably aligned in the a-b plane 
and neglect terms proportional to J 4 in comparison to 
those with J 2 . The constant a must be positive to fix J 
in the a-b plane. 

The fourth-order terms in Eq. (|l|) lift the 0(2) degen- 
eracy in favor of the four-fold symmetry. The Hamilto- 
nian (Q) can be treated quasi-classically since J x and J y 
are rather large (J=8, 15/2, and 6 for Ho, Er, and Tm, re- 
spectively). We introduce the angle (f> for the orientation 
of the moment J in the a-b plane, so that J x — J cos(f> 



H 



~\W + bJ\l + cos(40)) - /icos(0 - fa), (2) 



where h is the absolute value of the external magnetic 
field multiplied by 5 J/ig /4, fa is the angle between direc- 
tion of magnetic field and reference tetragonal axis {e.g. a 
axis). Exploiting again a large value of J, it is natural to 
assume that 2bJ 4 ^§> a. Then the potential energy in (Q) 
has deep wells at directions = ±7r/4, ±37r/4. Non- 
zero magnetic field destroys the tetragonal symmetry, but 
the symmetry breaking is small unless h exceeds the value 
6 J 4 corresponding to fields essentially large than 20 T [Jt) . 
Thus, initially continuous moment J is reduced to a dis- 
crete variable taking only four values. This is a kind of 
so-called clock model with 4 positions of the "hand" . 

Let us denote 4 states of the "hand" as where 
k = 0, . . . , 3. Then the Hamiltonian (|J) is equivalent to 
4x4 matrix 



~Hkk' = £k8kk> + U>(Sk,k'-l + 5k,k'+l)i 



(3) 



where ty. = — /icos((2fc + 1)tt/4 — 4>h) and if is a matrix 
element of the "kinetic energy" between adjacent 

states of the "hand" . In the approximation h <§C bJ A all 
the overlap integrals are the same (w). The diagonal- 
ization of the matrix (j^) can be performed explicitly at 
any values of the parameters h, fa, and w. The lowest 
eigenvalue gives the energy of the ground state 



E : 



\j \ (lw 2 + h 2 + v/(4w 2 + h 2 ) 2 - h 4 cos(2(/., 1 ) 2 ) ■ 

(4) 

The result (Q) explains the orientation dependent satu- 
ration in high magnetic field |3|,|| (see Fig. Q). Namely, 
it has been observed in M that the magnetization along 
magnetic field reaches saturation values depending on di- 
rection according to the law M s (fa) cx cos(7r/4 — fa). 
The moments in such a state are directed along an 
easy direction closest to the direction of the field. The 
direction-dependent saturation Fig. ^ is an intermediate 
asymptotic corresponding to the field interval w <C h <C 
bJ 4 . At higher fields h ^> b,J A the direction independent 
saturation is reached. 
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Further we employ this simplified description of the 
localized moments with minor modification to incorpo- 
rate the collective effects. Let us introduce variables 
k r = 0, . . . , 3 at each site of the lattice and the inter- 
action Hamiltonian 



Tii =^~~^ K r r i cos — (k r — k r r) i r ,r' cos7r(/c r — k r i), 



(5) 

where K T y and L r y are interaction energies. Only near- 
est neighbor interaction will be taken into account, so 
that the only non-zero constants are if r r ± a = if r , r ±b = 

— f-K"||) ir,r±a = £r,r±b = _ 2^ll' ^r,r±c/2 = K± and 

L T>T ± c /2 = Lj_, where a, b,c are the primitive vectors of 
the lattice. We assume that ify > K±_ > 0, \L\\\ > \L±\. 
For simplification we put (L\t = L± = 0). 

In this article we find only the ground state of the total 
Hamiltonian Ti = Ti s + Tii where the single-ion Hamil- 
tonian TL S is the direct sum of matrices (||). The ground 
state will be supposed to be homogeneous in-plane. The 
existing experimental data have no agreement upon the 
in plane structure of HoNi 2 B 2 C. The experiments per- 
formed on single crystals by Goldman, at al. have 
clear presence of the satellites at 0.585a* in the tempera- 
ture range of 5-6 K. On the other hand neutron powder 
diffraction measurements by Grigereit, at al. did 
not find evidence of these satellites. Both groups agree 
that the structure in the a-b plane transits to the ferro- 
magnetically aligned sheets below 5 K. We will assume 
ferromagnetic ordering in the a-b plane. This implies 
that Km > 0. The antiferromagnetic order in c-direction 
implies that K± > 0. 

We apply a variational procedure, assuming the total 
wave function to be a product of single-ion wave func- 
tions. The wave functions at sites belonging to the same 
a-b plane are assumed to be identical. Neighboring planes 
are assumed to belong to different sublattices 1 and 2. 
The single-ion states are described by two set of varia- 
tional parameters a.k,(3k, k = 0, ... ,3 where a>k and (3k 
are the amplitudes for a "hand" to be in the fc'th poten- 
tial well for the first and second sublattice respectively. 
With this premises the energy per site is 



TL V 
if 



,„ - ^[(ag - alM - ft) + (a\ - aj)(ft - ft)} 



□2\2l 



[H - K - «S) a + (ft - m 2 + - Pi) 



+w[(a + a 2 )(ai + a 3 ) + (A> + fo)(Pi + AO] 
-y[(«§-«2) + ($-#)] 
- al) + ((31 - ft)}. 



(6) 



where h x — hcos(ir/4 — (p^) and h y — hsm^/A — (j)^), 
due to the symmetry the consideration is constricted to 



the range < <f>h < tt/4. The minimization of the Hamil- 
tonian (H) with the constraints E a ? = Eft 2 = 1 leads 
to a system of nonlinear equations which allows a partial 
separation of variables. Proper variables are determined 
by following relations 

cos Ci = chq — a\ — «2 + a\, 
cosQ 2 = 2 o -ft-ft + ft, 

2,2 2 9 

cos rji = a + a 1 — a 2 



cos m = ft + (3\ - ft 



"35 
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In terms of Q and rji the minimization conditions read as 
follows 





_ cos Ci 
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— if 11 cos 771 


— 2 cot 771 


= h x 


+ hy, 



(8) 



where w has been set to 1. The systems of equations (g) 
have been solved numerically 

In the general case < <j>h < 7r/4 there are three phases 
and two first order phase transitions at fields h c i((f>h) and 
h C 2(<ph)- I n the range of small magnetic field the most en- 
ergy favorable is antiferromagnetic phase ( AF) , in which 
spins in one sublattice are parallel whereas spins of the 
neighboring sublattices are antiparallel. At increasing 
field it is replaced by the ferrimagnetic phase (FM), in 
which spins of the neighboring sublattices are perpen- 
dicular. Finally the spin-flip proceeds into a phase with 
almost parallel spins in the neighboring sublattices (para- 
magnetic or spin-flip phase (PM)). In special cases of 
4>h = and 4>h — the numbers of phases and phase 
transitions are reduced by one. Namely, for <ph = the 
system undergoes a transition from the AF into the FM, 
final paramagnetic saturation happens for experimentally 
unreachable magnetic field h > bJ 4 . For <ph = 7r/4 the 
system takes up from the AF directly into the PM. Nu- 
merically found diagram in h-tfih plane is presented on 
Fig. I for K± = 3, if || = 4. 

For the two phase transitions and magnetization dis- 
continuities at the transition points we find 



h c i 



cos <p h 



7-c2 



sin <ph 



AMi = M c cos < 



AM 2 = M c sin <j> h 



(9) 



(10) 



where h c = h c (Kx,K\\) and M c = M c (Kj_, ify). The 
graphs of h c and M c vs. Ku for some particular values of 
if_L are given at Figs. || and |^, respectively. The graphs 
of magnetization vs. magnetic field for different direction 
of the field are shown in Fig. ||. 

A series of simple relationships can be obtained in a 
limiting case of ify 3> w, well approaching the condi- 
tions of the experiments. Then the "hands" are tightly 
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bound to the tetragonal easy directions. Magnetic field 
lifts four-fold symmetry and one of the "hands" (labeled 
by oti for definiteness) is set to the easy direction closest 
to the direction of magnetic field ccq = 1 and on = 0, 
i = 1,2,3. The other "hand" can take three positions 
(03 = 0) depending upon strength and direction of mag- 
netic field. The three positions correspond to the three 
phases, i.e., AF 0\ = 1, # = 0, i = 0,1, FM 0\ = 1, 
0i = 0, i = 0,2, and PM 0$ = 1, 4 = 0, i = 1,2. The 
energies, magnetization and transition lines can be found 
straightforwardly resulting in 



h c (K ± ,K {l -oo) = ^, M C (K ±) K {] 



°° ) = 7r 

(11) 



These facts are confirmed clearly by numerical analysis 
(see Figs. || and |]). 

The results fit reasonably well the experimental facts. 
P. C. Canfield et. al. Q have observed three phase tran- 
sitions instead of two predicted by our theory. We denote 
the experimental lines as -ff c i(0/i), H C 2((j)h), and H c ^{<j)h). 
It has been found M that the dependencies are 



H c i — 



3.36kG 



cos(vr/4 - <j) h ) ' 



H C 3 = 



9kG 

sin 4> h 



There has been no evaluation given for H C 2- Our critical 
lines h c i and h c i plausibly correspond to the experimen- 
tal lines H C 2 and H C 3 respectively. There is a perfect fit 
between h c2 (see Eq. ([)])) and H c3 giving h c s» 9kG. We 
have already noted the excellent agreement between the 
theory and experiment for the angular dependence of the 
saturation magnetization. 

The drawback of the model is that one extra phase, ob- 
served in the experiment between lines H c \ and H C 2, is 
absent in our model. We presume that this is the chiral 
phase seen in the zero-field neutron diffraction experi- 
ments HQ (the satellite 0.915c*). Of course the spiral 
period might change in the presence of magnetic field. 
The spiral phase cannot be principally obtained in the 
framework of the two sublattice model. The number of 
layers should be at least equal to the ratio of the spiral 
structure period to the spacing between nearest holmium 
layers. The physical reason for appearance of the long- 
periodic modulated phase can be either RKKY forces 
or quantum tunneling (w) which becomes active when 
the main collective interaction K± favoring antiferromag- 
netism is suppressed by the magnetic field. Independent 
estimates of the quantum parameter w ||] show that it 
is rather small (w ~2-2.5 K). Nevertheless, it cannot be 
neglected since it determines the magnetic susceptibility 
in the low-temperature limit. 

We have shown that a simple 4-position clock model 
stemming from the quasiclassical single-ion Hamiltonian 
explains the angular dependencies of the saturation mag- 
netization and at least one of the transition points ob- 



served in the experiment. We believe that this fact indi- 
cates convincingly that the 4-position model is valid for 
the magnetic subsystem of RE nickel boride-carbides. 

Incorporating the simplest inter-spin interaction we 
have found two phase transitions between AF, FM and 
spin-flip phases. We expect that in a slightly modified 
model modulated phases are the ground-states in some 
range of parameters. 

Among the predictions which can be checked experi- 
mentally we emphasize following ones: 

i) In FM phase the magnetization vectors of the sublat- 
tices are perpendicular each other. 

ii) The jumps of magnetization are simple functions of 
direction (sec Eq. (10)). 

Our next purpose is the thermodynamics of the model 
as well as consideration of the long period chiral struc- 
tures. 

We are grateful to P. C. Canfield and collaborators for 
the opportunity to read their articles |^,|3| prior publi- 
cation and especially to D. Naugle for numerous discus- 
sions. 
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FIG. 1. Magnetization vs. magnetic field for magnetic 
field direction (from top to bottom) (f>h = 45, 36, 27, 18, 9, 
and degrees, respectively, the single-ion model (see 
Eqs. ([|^])). 4>h — corresponds to the direction [100] and 
h = 7r/4 to [110]. 



FIG. 2. h~4>h phase diagram, Kj_ 
sublattice model. 



3, K\\ = 4, the two 
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FIG. 3. h c (see Eq.(g)) vs. K\\, for interplane interaction 
(from top to bottom) K± = 3, 1.5, 0.5, and 0.2, respectively, 
the two sublattice model. 

FIG. 4. M c (see Eq.((l|)) vs. K ]} , the two sublattice 
model. 

FIG. 5. Magnetization vs. magnetic field, K± — 3, 
K\\ = 4, the two sublattice model. 
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Fig.3 



